arXiv:1509.02605v2 [math.DS] 15 Apr 2016 


Collision Index and Stability of Elliptic Relative Equilibria in 

Planar n-body Problem 

Xijun Hu^* Yuwei Ou^ ^ 


^ Department of Mathematics, Shandong University 
Jinan, Shandong 250100, The People’s Republic of China 
^ Chern Institute of Mathematics, Nankai University 
Tianjin 300071, The People’s Republic of China 


Abstract 


It is well known that a planar central configuration of the «-body problem gives rise to solutions 
where each particle moves on a specific Keplerian orbit while the totality of the particles move on a 
homographic motion. When the eccentricity e of the Keplerian orbit belongs in [0,1), following Meyer 
and Schmidt, we call such solutions elliptic relative equilibria (shortly, ERE). In order to study the linear 
stability of ERE in the near-collision case, namely when 1 - e is small enough, we introduce the collision 
index for planar central configurations. The collision index is a Maslov-type index for heteroclinic orbits 
and orbits parametrised by half-lines that, according to the Definition given by authors in m, we shall 
refer to as half-clinic orbits and whose Definition in this context, is essentially based on a blow up 
technique in the case e = 1. We get the fundamental properties of collision index and approximation 
theorems. As applications, we give some new hyperbolic criteria and prove that, generically, the ERE of 
minimal central configurations are hyperbolic in the near-collision case, and we give detailed analysis of 
Euler collinear orbits in the near-collision case. 
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1 Introduction 


For n particles of mass mi ,..., let qi ,..., q„ € be the position vectors, pi, ...,pn e be the momentum 
vectors. Setting dtj =|| qi - qj ||, the Hamiltonian function has the form 



( 1 . 1 ) 


U will be defined on configuration space 
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where A = {x € 1.^” : 3/ t j,Xi = xj] is the collision set. A central configuration is a solution (^i, ...,qn) - 
of 


- Arrijqj - 



qn) 


( 1 . 2 ) 


for some constant A. An easy computation shows that A = U{a)II{a) > 0 where I{a) = ^ mj\\aj\\^ is the 
moment of inertia. Otherwise stated, a central configuration with /(a) = 1 is a critical point of the function 
U constrained to the set £ = {x ^ K \ I{x) - 1). 

It is well known that a planar central configuration of the n-body problem gives rise to solutions where 
each particle moves on a specific Keplerian orbif while fhe fofalify of fhe particles move on a homographic 
motion. Following Meyer and Schmidl lISTl . we call these solutions as elliptic relative equilibria and in 
shorthand notation, simply ERE. Specifically, when fhe eccenfricify e = 0, fhe Keplerian elliptic motion 
becomes circular and all the bodies move around the center of masses along circular orbits with the same 
frequency. Traditionally these orbits are called relative equilibria . 

As pointed out in OTl . there are two four-dimensional invariant symplectic subspaces, E\ and E 2 , and 
they are associated to the translation symmetry, dilation and rotation symmetry of the system. In other 
words, there is a symplectic coordinate system in which the linearized system of the planar n-body problem 
decouples into three subsystems on E\, E 2 and £3 = (£1 U£ 2 )'‘', where X denotes the symplectic orthogonal 
complement. A symplectic matrix AI is called spectrally stable if all eigenvalues of AI belongs to the unit 
circle U of the complex plan, while AI is called hyperbolic if no eigenvalues of AI are on U. The ERE is 
called hyperbolic (resp. stable) if the monodromy matrix AI is restricted to £3, that is, AIl^j is hyperbolic 
(resp. stable). 

There are many interesting results for the linear stability of ERE (cfr. |[^ [34l |38j [33 and references 
therein). Many of them, investigated the relative equilibria for e small enough and as far as we know, only 
few of them studied the linear stability of ERE with e € [0,1). To our knowledge, the elliptic Eagrangian 
solution is the only case that is well studied. The Eagrangian solution which was discovered by Eagrange in 
1772 It22l is the ERE of the equilateral triangle central configuration in the planar three body problem. 

It is well known that the stability of elliptic Eagrangian solutions depend on the eccentricity e and on 


P = 


21{m\m2 + mitriT, + m2m^) 
(nil +1712 + 


(1.3) 


Eong et al. used Maslov-type index and operator theory to study the stability problem, and gave out a full 
describe of the bifurcation graph (cfr. ifTTll . lfTSl ). Moreover, Wang et al. built up a trace formula for linear 
Hamiltonian systems and Sturm-Eiouville systems, and used it to give an estimate of the stability region as 
well as of the hyperbolic region I^ . EOll . 

In the study of the near-collision case, that is, when 1 - e is small enough, a blow-up technique from R. 
Martinez, A. Sama, C. Simo Il28ll is very useful to carry over our analysis. The authors considered 4D linear 
system depending on a small parameter cr > 0 and the singular limit for cr ^ 0 . Based on it, they computed 
the trace tr\ = Ti - 1 - / 1 “\ tr 2 = A 2 + where Ai,Aj^, i = 1,2 are the eigenvalues of monodromy matrix. 
Under a “nondegenerate condition” they describe the asymptotic behaviour of log |tr,|, i - 1,2 and tr 2 . Their 
study includes the ERE of Eagrangian equilateral triangle and Euler collinear central configurations. 

Motivated by the these results, we will use blow-up technique and index theory to study the stability 
problem of ERE. The index theory we used will be based on the Maslov-type index. The Maslov index is 
associated to a given continuous path of Eagrangian subspaces and the Maslov-type index is assigned to a 
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path of symplectic matrices. We briefly review the Maslov index theory in §2.2 and give its relation with 
the Maslov-type index. For reader’s convenience, we now roughly describe the Maslov-type index theory. 
Let Sp(2n) be the set of symplectic matrix in equipped with the standard symplectic structure, and set 
hn be the identity matrix on . Let y e C([0, T], Sp(2n)) with 7(0) = hn, for m € U, roughly speaking, 
the Maslov-type index icj{y) is the intersection number (by a small perturbation) of y and {Me 

Sp(2n),det(M - ojhn) - 0). Please refer to 1241 for the details. 

By the blow-up technique, the limit of ERE, as e ^ 1, can be described by two heteroclinic orbits Zq, I+ 
connected P+. (Cfr. to Figure [^. Throughout of the paper, we denote by y^ the fundamental solution of the 
essential part of ERE, that is, ye(t) = JS{t)ye(t),t e [0,2;:], 7^(0) = id, where S{t) is defined in Equation 
( |2.2| ). Given a symmetric matrix R we shall denote by /li(R), the smallest eigenvalue of R. When the limit 
equilibrium P± is non-hyperbolic, we have the following result. 


Theorem 1.1, Let ao ^ S be a planar central configuration which satisfies 


AfiD^Ulsiao)) < -luiao). 


(1.4) 


where D^Ulsiao) is the Hessian of U restricted to S at ao. Then Zi(ye) oo as e ^ 1. 

Eor p^AiiD^Ulsiao)) = ~ n, let s - j min {, 1/8), e - Ife < e^, then we have 

Vn /eS 

Zi(y.)>2^1n(— -6. 


(1.5) 


Erom Il241 . it is well-known that, for any m e U, \icj - Zi| < n. Then Theorem 1.1. shows also that 
i(xji7e) —> oo as c —> 1, which implies there exists a sequence efico) converging to 1, such that the system is 
m-degenerate. 

It is well known that any T -periodic solution is a critical point of the action functional 


r{q)^ 


r*l '' 

fJW-" 


wmf 


U{q)^dt 


defined on loop space W^’^(R/rZ, A). Eet now Xe be the ERE corresponding to ao with eccentricity e, and 
let (pixe) be the Morse index of Xe (meaning that it is the total number of the negative eigenvalues of T'"{xe)). 
Since the Morse index is equal to Maslov-type index (cfr. to Eemma [53| We have f{Xe) is not less than the 
Maslov-type index Zi(ye) coming from the essential part. Theorem | 1.1 [ implies that, if ao satisfied ( |L4| ), then 
(piXe) —> oo as c —> 1 m. 

The above theorem is related to the result of the interesting paper of V. Barutello and S. Secchi 
They defined a collision Morse index for one-collision solution in n-body problem with a homogeneous 


potentials, and proved that the collision index is infinite under the condition (1.41 for the Newton potential. 
Their results show that a one-collision solution asymptotic to ao which satisfied ( |L4| ) cannot be locally 
minimal for the action function. A Morse-type index theorem both for colliding and parabolic motions, will 
be given in ['5lj. 

The most interesting case, however, is precisely when P± is hyperbolic, in this case we can define the 
Maslov index for heteroclinic orbits and half-clinic orbits. (Cfr. Equations ( |3.2| ) and ( |3.4| )). Eor half-clinic 
orbits, we mean a solution of Hamiltonian system x{t) defined on R"'' or R“, where R'*' and R“ stands for 
the non-negative and non-positive half-line, respectively and such that the initial condition x(0) belongs to a 
Eagrangian subspace whilst x{t) converge to an equilibrum point when t —> ±oo. Also in this last case , we 
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assign a Maslov index to both /q and 1+ and we shall refer to as collision index. After defined the collision 


index, we shall prove Theorem 3.3 and we shall refer to as approximation theorem. Let us now show that, 
under a suitable non degenerate conditions for e ^ 1, the Maslov index for is convergent to the sum of 
collision index on Zq, l+. This is a main part. (Cfr. §3.1 for the details). In the study of stability problem, the 
Dirichlet, Neumann, periodic, anti-periodic boundary condition play an important role. Our key idea is to 
use the Maslov index corresponding to these 4 kinds of boundary conditions for determining the stability. 

Throughout of the paper, we always let Lj = R-^ 00, Vi = O0R-^ be the Lagrangian subspace in (R^^ a>o) 
which corresponding to the Dirichlet and Neumann boundary conditions respectively, and we always omit 
the subscript if no confusion is possible. For ERE, by using the approximation theorem, we get the following 
result. 

Theorem 1.2. Let ao e S be such that /li(D^Z7|£(ao)) > -|Z7(ao) and we assume that ao is nondegenerate 
and collision nondegenerate. Let (pla^) be the Morse index of a^ which is the total number of negative 
eigenvalues of D^U\£,{a{y). For \ — e small enough, we have, ye{2n)Vd iti Vd, 


b(yd, 7 e{t)Vd, t € [0, In]) = k + i{Vd\l+). 


(1.6) 


and ye{2n)V„ iti V„, 


bi^n, 7e{t)Vn, t € [0, 27r]) = 2(p{ao) + i{Vd\ z+), 


(1.7) 


where itl means transversal, k = 2n - 4, i(Vd', 1+) is the collision index on 1+ defined by {3.4) and // is the 
Maslov index. 


The definition of collision nondegenerate index is given in Definition |3.2[ The degenerate problem on Zq 
will be discussed in §3.2. We observe that, in contrast with respect to the nondegeneracy condition along Zq, 
we didn’t establish a useful criterion for detecting the nondegeneracy along Z+. 

If the central conhgurations have brake symmetry (cfr. Dehnition |4.1[ ), the collision index of heteroclinic 
could be decomposed into the sum of index on half-clinic orbits and this will simplify the computation. 
To our knowledge, the Eagrangian and Euler central conhgurations both have brake symmetry. Another 
example is the l+n central conhgurations, that is regular polygon conhgurations with a central mass. It will 
be interesting to provide central conhgurations without this symmetry property. 

As an application, we study the stability of ERE for minimizer central conhgurations. Eor a central 
conhguration ao, it is obvious that D^U\siao) has a trivial eigenvalue 0 which comes from the rotation 
invariant. The central conhguration ao is called nondegenerate minimizer if all the nontrivial eigenvalues 
are bigger than 0, while ao is called strong minimizer if all the nontrivial eigenvalues are bigger than U{ao). 


Theorem 1.3. 'We assume that ao is a nondegenerate minimizer that satisfies the collision nondegenerate 
condition. lf\ — e is sufficiently small, then the ERE is hyperbolic. 


In the case e = 0, Moeckel conjectured f2\ that a relative equilibrium is linearly stable only if it as¬ 
sociated to a minimizing central conhguration. Our results show that in the case that \ - e small enough, 
it is generally hyperbolic. By the way, we conjecture that Theorem 1.3 is true also without the collision 
nondegenerate condition. 

In the case ao strong minimizer the following result holds. 


Theorem 1.4. The ERE of a strong minimizer ao is hyperbolic for any e € [0,1). 
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A typical example of nondegenerate minimizer central configurations is the Lagrangian central config- 

was 


urations, which is strong minimizer if yS > 8 . For Lagrangian orbits, the conclusion of Theorem 1.3 


proved in llT7]l without the collision nondegenerate condition and the result of Theorem 1.4 was proved by 
Ou [O. Another easy example is the 1 + 3-gon central configurations, that is, the regular triangular config¬ 
urations with a central mass. The three unit masses with unit distance away from the mass nic at the origin. 


As a direct application of Theorem 1.3 and|1.4[ we get the next result. 


Corollary 1.5. Let ao be a I + 3-gon central configuration with central mass me, for me e [0, ^^ 249 ^ ) 
and we assume that ao is collision nondegenerate. lf\ — e is sufficiently small, then the ERE is hyperbolic 
Eurthermore, if me € [0, then the ERE is hyperbolic for any e € [0,1). 

Another conjecture of Moeckel ||3 states that a relative equilibrium is linearly stable only if it has a 
dominant mass. For example the Lagrangian orbits and ERE of 1 - 1 - 3-gon relative to a strong minimizer 


have no dominant mass. Thus Theorem 1.4 can be considered as a support of Moeckel’s conjecture in the 


case of c > 0, so we guess Moeckel’s conjecture is also true in the case of ERE. 

The collision index plays an important role in the study of the stability problem. We shall give some 
conjectures for the collision index which are related with Y. Eong’s conjecture for the Maslov-type index of 
ERE. (Cfr. Remark [53] for further details). 

As a further application, we consider the ERE of Euler collinear central configurations iTnil . that we 
simply refer to as elliptic Euler orbits. The linear stability of this kind of orbits depends upon two parameters, 
e and 6, where the last one 6 € [0,7] only depends on mass mi, m 2 , m^. (Cfr. Appendix A of 1(281 and l(27l 
for further details). To our knowledge, the near collision case was firstly studied by R. Martinez, A. Sama, 
C. Simo |[28l . Y. Eong and Q. Zhou used Maslov-type index theory in order to describe the +1-degenerate 
curves; they also analysed the stability problem. It is worth noting that by their methods is not possible to 
explain the limit property of ±l-degenerate curves numerically proved by R. Martinez, A. Sama, C. Simo 
||2^ . Please refer to Eigure|7]and Eigurej^ Using the collision index, we explain the limit property. We show 
that 6 > 1/8 is equivalent to condition ( |L4| ). Theorem 1 1.1 1 implies the ±l-degenerate curves don’t intersect 
[1/8,7] X 1. In the case, S € (0,1/8), the collision index is well defined, we analysis the near-collision 
phenomena by the collision index. We can compute in detail for collision index on fi, but unfortunately, we 
can’t determine the collision index on /+ by analytical method. Instead, we develop a numerical method to 
compute the collision index. Based on numerical results A, the collision index strictly proved the behaviour 
of the ±1 in the near-collision case. Please refer to §5.2 for the details. 

This paper is organized as follows. We review the Meyer-Schmidt reduction and Martinez, Sama,Sim 6 
blow up technique at §2.1. We give a brief introduction of the Maslov index theory and we prove Theorem 
o in §2.2. The definition of collision index is stated and the approximation theorem is proved in §3.1. 
Some basic properties of the collision index are given in §3.2. In §4, we study the case of brake symmetric 
central configurations. The computation of the collision index along lo,l+ is given in §3.2, §3.3. We give 
some applications in §5. In §5.1, we study the minimizing central configurations and we prove Theorem 


1.3 and Theorem 1.4 We use the collision index to analyse the Euler orbits in §5.2. At last, for the reader’s 


convenience, we give the details of the numerical method used to compute collision index in § 6 . 


2 Blow up and limit index for the non-hyperbolic case 

This section includes some basic preliminaries. We first briefly review the decomposition of ERE by fol¬ 
lowing authors in OTII and the blow-up technique of Martinez, Sama and Simo 11291 in section §2.1, then we 
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review the fundamental property of Maslov index in §2.2, and give the proof of Theorem|l.l| 


2.1 Reduction and blow up method 

In 2005, Meyer and Schmidt strongly used the structure of the central configuration for the elliptic La- 
grangian orbits and symplectically decomposed the fundamental solution of the elliptic Lagrangian orbit 
into two parts, one of which corresponding to the Keplerian solution and the other is the essential part of 
the dynamics, needed for studying the stability. For the reader’s convenience, we briefly review the central 
configuration coordinates, by following Meyer and Schmidt lISTTl . 

Suppose that Q = {qi ,..., qn) € with mass nii, ...,m„ is a central configuration, and P = {pi, ...,Pn) e 

R^”. Let Ij be the identity matrix on W, J 2 j = [ |. We denote by J„ = diag{J 2 ,J 2 ) 2 nx 2 n and 

V h Oy ) 

M - diag{mi,m\,m 2 ,m 2 , ...,m„,m„) 2 nx 2 n- We assume that t i-> x(0 is a periodic solution of ERE, then the 
corresponding fundamental solution is 


7(0 - J4nH"{x{t))y{t), 7(0) = hn. 


( 2 . 1 ) 


As in lISTl Corollary 2.1, pag.266]. Equation ( |2.1[ ) can be decomposed into 3 subsystems on E\, E 2 and 
£3 = (El U E 2 )'^ respectively. The basis of Ei is (0, m), {Mu,0), (0, v), (Mv,0), where u = (1,0,1,0,...), 
V = (0,1,0,1,.), andE 2 is spanned by (0,<3), (M<a,0), (0, J„<3), (J„M<a,0). EorX - (g,z,w) € R2xR2xR 2«-4 
and Y = {G, Z, W) e x R^ x we consider the linear symplectic transformation of the form Q = 

PX,P - P-^Y, where P is such that JP = PI,P^MP - hn (lEIl, p263). Now B{t) - H'\x{t)) in this new 
coordinate system has the form B{Q) - Bi® B 2 ® P 3 , where B, = The essential part Bi,{t) is a path of 
(4n - 8) X (4n - 8) symmetric matrix. 

In the rotating coordinate system and by using the true anomaly as the variable, Meyer and Schmidt |[^ 
gave a very useful form of the essential part 


m) = 


h -4/2 

T T — I ’ 

^k/2 h l-Hpcos(?) / 


t € [0,27r], 


where k = 2n - 4 and e is the eccentricity, t is the true anomaly and 


D = jP'^D^U{Q)P\ 


wsR*’ 


with A - 


ujQ) 

HQ)' 


( 2 . 2 ) 


(2.3) 


We denote by P := 4 + 2), which can be considered as the regularized Hessian of the central configurations. 
In fact, for ao € S which is a central configurations, then I{ao) = 1. With respect to the mass matrix M 
inner product, the Hessian of the restriction of the potential to the inertia ellipsoid, is given by 


D^U\s{ao) = M~^D^U{ao) + U{ao). 


Then we have 


p-^D^U\e{ao)P = P^D^U{ao)P + U{ao), 


(2.4) 


and thus 


R = 


1 

f/(ao) 




(2.5) 
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Thus the corresponding Sturm-Liouville system is 


-y- IJk/iy + 


R 

\ + e cos(t) 


y-0. 


( 2 . 6 ) 


In order to study the singular limit case e ^ 1, we use a change of coordinates as in Il29l . In fact, in the 
case of Newtonian potential, this can be interpreted as a McGehee change of coordinates (cfr. Il30ll . Il32]| ). 
Let q = (l+e cos(t))^^^, Q = -2q and change the time variable to t, where dt = qdr. Throughout the paper, 
we always use x' = ^ and x = Then we have 


q' = -2^Q’ 


Q' = 2Q^ + q^ 


1 . 


(2.7) 


We observe that (2.7 1 is well defined for ^ = 0 and its first integral is £ = q + y ~ 1)- An easy 
computation shows that, for the orbits with eccentricity e, the first integral with E = = -e. The system 

has two equilibria P+ = (0, + V2) lying on the level set E = 0. We distinguish the level set £ = 0 into two 
orbits 


Zo = {(^,0€]r2|^^O,|(2I< V2}, (2.8) 

and 

l^ = {{q,Q)^^\>0,Q^+q^ = 2). (2.9) 

On Iq, we have 

qioiT) = 0, 00(t) ^ - V2tanh(^T), (2.10) 

and the system on Z+ is 

ql=-\qQ, Ql=-^- (2-11) 

The solution is given by 

^/^(r) = V2/cosh(^^), 0 _^(t) = V2tanh(^^). (2.12) 

For convenience, we also let I~ = {{p, Q) ^ Iq,Q < 0) and = {{p, Q) e h,Q > 0), similarly, let 
K = {(P, Q) e Z+, 2 < 0) and Z+ ^ {{p, Q) €l+,Q> 0}. Obviously, Zq ^ U Zg and Z+ ^ Z; U Z+. 

Throughout of the paper, we always let ye (e € [0,1)) be the fundamental solution of ( |2.2| ). For simplic¬ 
ity, je can also be considered as a function of r. For q ^ 0 we consider the matrix S = diag{q^Ik, q^^h) e 
Sp(2Z:) and for T = t{2ti), we have S{T) = 5(0). We let jeix) = 5(T)ye(T)5“^(0) and we observe that the 
associated monodromy matrix is similar to the one of y,,. A direct computation shows that 

^ye = JBye, ye{0) = hk, r€[0,r], (2.13) 

dT 

with 

^ _ E ~ q2k/2 

jlk + q^ki2 q^h-R ] 

1 


(2.14) 









Figure 1: Phase portrait of (2.7 1 from 


Figure 2: An illustration of the sections used in the 


proof of Theorem 3.3 


The linear system ( |2.13 1 is well-defined also when e = 1. In this case, F = 0, the system has two 
equilibria corresponding to P+, and the system can be considered as two heteroclinic orbits. 

Proposition 2.1. P+ is hyperbolic if and only if AfR) > -f- 


Proof We observe that at points P+, the linear part with form D± - JA 


±^4 


±^4 

-R 


. P+ is hyperbolic 


if and only if the eigenvalue of D+ is not on the imaginary line. Since R is diagonalizable by choose suitable 
bases, the results is from simple computations. □ 


Given e < 1/8, we define fhe following sections (See figure 2) 

2^0 ^ {(<?, G)|0 <^<£,2 = 0), Xi = {{q, 2)|0 <q<s,Q^-^ + s], 

1^2 = {{q,Q)\q = e,--j2-E^ <Q< s - V 2 }, X 3 -((-?, 2)10 < V 2 -^<£, 2 - 0 ). 

2^4 ^ {{q, Q)\q ^e, V 2 -e< 2 < V2 - eA, S5 ^ {{q, 2)|0 <q<s,Q= V2 - e). 

If e = E(e) {{q, Q), E = -e] denofes fhe energy level sef. A direcf compulation show lhal E{e) 
inlersecfs E/, i = 0, • • • ,5 simply, i.e. infersecf af exacfly one poinf when e < s < 1/8. In fhis case, 
fhe Poincare map P, : E,_i i-> E,, / = 1, - • • ,5 is well defined. In facl, e < ensures fhe intersection 
wifh Eo,E3. For fhe intersection of E{e) and Ei, we observe lhal, since Q - - V2 - 1 - e, are solutions of 
q^{{- ^ + e)^ ll-\-q^ jl- 1) = -e, we gel q\ = \{e{1^I1- e) ± ((e(2 V2-e))^ - 8 e) 2 ). For e < 1/8, e < e?, 
we have q^_ < e^ < q\, which guaranlee fhe simple inlerseclion. Since E{e) is convex. Ibis also ensure fhe 
simple inlerseclion wifh E 2 , and we have fhe simple intersection of E 4 , E 5 by symmelry. 
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Following Il 29 l . we denote by Tig > 0 the time defined by Qigijig) = - ^|2 + s and r/^ > 0 such that 
qi^{-Ti^) - £. It is obvious Qig{-Tig) - - s and = s. Tig and are finite and independent 

of e once e is fixed. Let qo = q{ 0 ) and t\, T2 be the smallest positive time such that {q{Ti), Q{t\)) € Xi 
and {q{T2), Q(t2)) ^ '^2 ■ It is clear that qo, ti and T2 depend on e. Moreover qo ^ 0 , ti —> Tig and 
T jl - T2 ^ when e —> 0 . Similarly, let T4, T5 be the smallest positive time such that {q{T^), Q{t/\)) e S4 
and {q{T^), Q{ts)) € E5. We have T4 - T '/2 —> ti^ and T - T5 Tig. The next lemma was proved in |[ 29 ll : 
however, for the reader’s convenience, we proved it with a slight sharp estimates. 

Lemma 2.2. For s < \l%, e < e^, we have 
(a) 


(b) 


(c) 


V21n( 




) < T2 - Tl < 


ln(- 


V2-e ^(^1) 






\Q{t) + V 2 |(fT < 2e. 


(2.15) 


(2.16) 


(2.17) 


Proof. To prove (a), we observe that for r € [ti, T 2 ], - V 2 < 0 ( t ) < 


V2 

2 

rri 


get 5 ( V2 - £)q(T) < q'(T) < 
note that q{T)dT < q'{T)dT < 

To prove (c), let F{q) - v /2 - ^J 2 
have 


V 2 + £. Multiplying by -^(t)/2 we 

q{T) and hence Equation (|2.15|) by a direct integration. To prove (b), please 
2 


V 2 - 


-q{T2). Equation (|2.16|l follows by observing that q{T2) - s. 

Erom (b), we 


q^. Then F{q) < cq for q < e, with c - 


Ice 


V2+ V2-£2 ■ 


f"F(q)dT< . 

Jti y2- e 


(2.18) 


Please observe that, near P_, 1+ is graph of F{q) - V2. Eet y - Q + V2 - F{q), then 0 < y < £. By a direct 
computation it follows that 


y-( 


2V^ 


^J2- q^)y + y^/2 = - V2y(l + oi). 


where |oi| < e. Thus we have - V2y(l + £) < y' < - V2y(l - £), then 


r 

Oti 


1 £ 
ydT < — - y(Ti) < 


Erom ( |2.18||2.19| ), we have 

-T2 


r ' le + v2i^/t < r 

Jti Jti 


V2(1-£) V2(1-£) 


£ 2,C£ 

(y + F(q))dT < —- 1 -- < 2e. 


V2(l - £) V2 - £ 


(2.19) 


( 2 . 20 ) 

□ 
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2.2 The index limit on the non-hyperholic case 

We start by briefly reviewing the Maslov index theory [31[33| in this subsection. Let be the 

standard symplectic space, and Lag{2n) the Lagrangian Grassmanian, i.e. the set of Lagrangian subspaces 
of For two continuous paths L\{t),L 2 {t), t € \a,b\ in Lag(2n), the Maslov index /r(Li,L 2 ) is an 

integer invariant. Here we use the definition from |(6l. We list several properties of the Maslov index. The 
details could be found in |(6l. 

Property I (Reparametrization invariance) Let g : [c, d] \a, b\ be a continuous and piecewise smooth 
function with q{c) = a, g{d) = b, then 

li{U{t),L 2 {t)) - /r(Li(p(T)),L 2 (p(T))). (2.21) 

Property II (Homotopy invariant with end points) For two continuous family of Lagrangian path L\{s, t), 
L 2 {s, t),0 < s < \, a < t < b, and satisfies dimL\{s, a) D L 2 (s, a) and dimLi{s, b) D L 2 {s, b) is constant, then 

^i{U{0,t),L2{0,t)) = gi{Li{\,t),L 2 {l,t)). (2.22) 

Property III (Path additivity) a < c < b, then then 

L2(t)) = L2{t)\[a,c]) + L2{t)\[c,b])- (2.23) 

Property IV (Symplectic invariance) Let y(l), t € \a, b^ is a continuous path in Sp(2n), then 

L2(0) - gi{y{t)Li{t), y{t)L 2 {t)). (2.24) 

Property V (Symplectic additivity) Let Wi, i = 1,2 be symplectic space, L\,L 2 € C{{a,b^,Lag{Wi)) and 
Li,L 2 £ C{\a,b^,LagiyV 2 )), then 

/r(Li(0 © Li(0, L2{t) © 1,2(0) ^ biLiit), L2(t)) + -^2(0)- (2-25) 

In the case L\{t) = Vq, L{t) - y(t)V, where y is a path of symplectic matrix we have a monotonicity 
property (cfr. Jill). 

Property VI (Monotone property) Suppose for j = 1,2, Lj{t) = yj(t)V, where yj{t) = JBj{t)yj{t) 
with yj{t) - hn- If B\{t) > 82 ( 1 ) in the sense that B\{t) - 82 ( 1 ) is non-negative matrix, then for any 
Vo, Vi € Lag{2n), we have 

M^o,ri^i)>M^o,r2Vi). (2.26) 

One efficient way to study the Maslov index is via crossing form introduced by 1371 . For simplicity and 
since it is enough for our purpose, we only review the case of the Maslov index for a path of Lagrangian 
subspace with respect to a fixed Lagrangian subspace. Let A(t) be a -curve of Lagrangian subspaces with 
A(0) = A, and let V be a fixed Lagrangian subspace which is transversal to A. For v e A and small t, define 
w{t) € V by V -I- w{t) e A(t). Then the form 

e(v)=4 0 j{v,w{t)) (2.27) 

dt t=o 

is independent of the choice of V (cfr. ll37l '). A crossing for A(t) is some t for which A(f) intersects W 
nontrivially, i.e. for which A(t) € Oi(W). The set of crossings is compact. At each crossing, the crossing 
form is defined to be 


F(A(0, W, t) - 2|A(/)nw- 


( 2 . 28 ) 
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A crossing is called regular if the crossing form is non-degenerate. If the path is given by A{t) - y{t)A with 
7(0 € Sp(2n) and A € Lag{2n), then the crossing form is equal to Jy{t)v, v), for v € 7(0“^(A(0 n W), 

where (, ) is the standard inner product on R^”. 

For A(0 and W as before, if the path has only regular crossings, following ll25l . the Maslov index is 
equal to 

^i{W, A(0) - m+(r(A(a), W, a)) + sign{r{A{t), IF, 0) - m-{r{A{b), W, b)), (2.29) 

a<t<b 


where the sum runs all over the crossings t e (a, b) and ni^, m~ are the dimensions of positive and negative 
definite subspaces, sign - ni^ - wT is the signature. We note that for a C'-path A{t) with fixed end poinfs, 
and we can make if only have regular crossings by a small perfurbafion. 

In confrasf with the definition given in equation ( 2.29| ), the Maslov index defined in OTl has the follow¬ 
ing form 

liRs{A{t), W) - ^sign{r{A{a), IF, a)) + ^/gn(r(A(0, W, 0) + ^sign{r{A{b), IF, b)). (2.30) 

a<t<b 


We observe that, for the non-degenerate path (i.e. L{t) n IF = 0 for f = a, b), 


Note that for M € Sp(2n), Gr{M) := {{x,Mx) \ x € R^”) is a Lagrangian subspace of the symplectic 
vector space (R^" 0 R^”, -m 0 to). Let y{t) be a path of symplectic matrices, A = Ai 0 A 2 e Lag{An), where 
A; € Lag{2n), for i = 1,2, then following ll37]l and by computing the crossing forms, we have 

/i(Ai 0 A2,Gr(7(0)) =MA2,7(0 Ai)- (2-31) 

For a continuous path y(t) € Sp(2n) with 7(0) = l2n, the Maslov-type index ic^iy) e Z is a very useful 
tool in studying the periodic orbits of Hamiltonian systems flM . The next lemma ( |[25l Corollary 2.1.) 
gives its relation with the Maslov index. 

Lemma 2.3. For any y{t), we have 

i\{y) + n= fi{A, Gr{y{t))), (2.32) 

and 

Liy) = loiGrito), Gr(ym, m € U\{ 11, (2.33) 

where A is the diagonal Gr{l 2 n), Gr{a)) = Gr{tol 2 n)- 

For Fi, F 2 € Lag{2n) and a Lagrangian path t 1 -^ A(t), the difference of fhe Maslov indexes wifh respecf 
fo fhe fwo Lagrangian subspaces is given in terms of fhe Hdrmander index, i.e. lITTlI (Th.3.5.) 

FFo,Fi;A(0),A(1)) -M^o,A)-MVi,A). (2.34) 


Obviously, 


FFo, Fi; A(0), A(1)) - Fi;c-'^''A(0),c-'^''A(1)), 


(2.35) 
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for e > 0 small enough. The Hormander index is independent of the choice of the path connecting A(0) and 
A(l). Under the non-degenerate condition, i.e. Vi, V 2 are transversal to A(0), A(l) correspondingly. Two 
basic properties are given below 

^(Uo, Vi;A(0),A(l)) = -^(Ui,l/o;A(0),A(l)), 


^(A(0),A(l);Uo, Vi) = -^(Uo,Ui;A(0),A(l)), 

If Vj = Gr{Aj), A{j) = Gr{Bj) for symmetry matrices Aj and Bj, then 

s{Vq, Ui; A(0), A(1)) = ]^sign{BQ - Ai) - 1 - ]^sign{Bi - Aq) - ]^sign{Bi - Ai) - ]^sign{BQ - Aq), (2.36) 

where for a symmetric matrix A, sign{A) is the signature of the symmetric form (A-, •). A direct corollary 
shows that 


|^(Uo,Vi;A(0),A(l))| <2n. 


(2.37) 


A sharp estimate for the difference of Neumann and Dirichlet boundary conditions has been given in |[26l . 

Let y be a fundamental solution of a periodic orbit, then y e C([0, T], Sp(2?i)) with y(0) = hn, as we 
have mentioned in the introduction, the Maslov index iJ.{Vn,yVn), )U(Urf,yU^) and Maslov-type index i\{y), 
i-\{y) play an important role in the study of stability problem. 

We come back to ERE, and we recall that ye{T) = S(T)y(,(T)S“^(0), with S = diag(^q2lk,q~^Ik) e 
Sp(2k). Please note that 5(7') = 5(0), and the path 5(t) is contractible in Sp(2n). In fact, if we set 
Sa = diag[q^lk, qa^h) with < 7 ^ = (1 -I- acos(t))^A^ then 5q; is homotopy to the constant path 5o(t) = hk by 
5Q. for a € [0, e\. We have 


Lemma 2.4. For e e [0,1), suppose dimGr{S a(7~)ye{T')S ^^{0)) n A is constant for a e [0,e], then 
p{Gr{ye), A) ^ p{Gr{ye), A). 


Erom Eemma|2.3|and Lemma 2.4 


we get 


Liye) = Life), Vm € U, ee [0, 1 ). 


(2.38) 


We first consider the Maslov index on L. Let 'To(t) be the fundamental solution on Iq, that is 


—T'o(t) = 7Bo'Po(t), T'o(O) = hk, T € [-CXD, +CX,), 
dr 


(2.39) 


with Bo - 


4 ^4 


r A 

V %4 -R 


Proposition 2.5. Suppose AfR) = -(1/8 - 1 - ri) with r\ > 0, then, we have 


p(Vd,^o(T)Vd,T e [0,To]) > [—to] 


(2.40) 


where [Z] denote the maximum integer which is not bigger than Z. 


12 





Proof. By changing the basis, we assume R = diag{A\, ...,Ak) where A\ < A 2 < ... < Ak and A\ < 
Based on the property V of the Maslov index, we have the decomposition 




(=1 


where 'B'(r) satisfies the equation 


^n(T) - JiBfP'oir), 4^'(0) = h, re [0, +cx,), 


(2.41) 


2/0 


with Bi = 


1 — 1 

Qi^ ^ ■ Since Bi\yi > 0, then r('PQ(T)V^, Vj, r) > 0, this implies that b(V^, 'Pp(T)l/j, r e 

V “ 4 “ ~Ai{R) ) 

[0, +To]) is nondecreasing with respect to tq. Moreover we have 

MV],'Fj,(T)V],Te[0,To])- y\Tj), 

0<Tj<To 

where v‘{Tj) = dim n '¥Q{Tj)V^. In order to compute the Maslov index //(Vj, r e [ 0 , tq]), we 

choose the basis ei = ( 1 , 0 )^ of and we let e[(T) = Then /i(V^, 'PQ(T)Tj, r € [ 0 , tq]) is equals to 

the number of zeros of /’{t) = det(M'(T)), for M\t) = (ci, Cj(t)). 
afr) bfr) 


Let 'P' (T) - 
equation 


C;(t) dfT ) 


. Then /'(r) = c,(t). From equation (2.41 1 , we get that c,(t) satisfies the 


V2 


^ I n tanh^(—r) - - + AfR) \ c,(t), 

Ci(0) - 0, c,(0) = 1. 

For i = 1, /li(/?) = - r\, then we have | tanh^(-^T) - 5 + /li(7?) < -r\. Using the Sturm comparison 

theorem, we know the number of zeros of c/(t) will be no less than This is complete the proof. □ 

In order to proof the Theorem o we need the lemma below. We will give an estimation of Maslov 
index on the period [ti,T 2 ]. Let e\ = min{ ^ , 1/8}, and y{T,T\) be the fundamental solution of (2.131 
with y(Ti,n) = hk, we have 


Lemma 2.6. For s < \s\, e < s^, 


B{Vd,y{T,T])Vd-,T € [ti,T 2 ]) > ^^ln(^)- 3 . 


(2.42) 


Proof. Let = 


k 




-R 


B-B- = 


|. Suppose I - q\ < E\ for some £1 < 1, then we have 
Ok 


h + q^kji 


^4 - ^4/2 ^ 


> -S\l2k- 


Let Bgj - B- - Eihk, by the monotonicity property of Maslov index, we have 

fi{Vd,y{T,Ti)Vd-,T e [ti,T2]) > fi{Vd,exp{{T - Ti)JBs)Vd-,T e [ti,T2]). 
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(2.43) 













Since R is diagonalizable, by using the o-product, we can split Bg, into the product of k two by two 
matrices, where the first factor (which is needed for computing the Maslov index) is given by Si 


1 - £1 


VI 

4 


i/8 + n-£i 


. By the direct sum property of Maslov index, we have 
yu(l/^,exp((T-Ti)7S£i)l/rf;T e [ti,T 2 ]) >/r(V],exp((T - ti)7Si)1/];t € [ti,T2]). 
Let/(£i,n) ^ sl-{r+ |)£i +n,S = diag{l, P ^ 



0 1 

f 1 f 1 

i 0 

(1 - £i)V 2 j 

[o 1 ) 


(2.44) 
, then 


We have 


P ' exp((T - t\)JB\)P - exp((T - t\)JB). 


MVrf,exp((T-Ti)7Sil/^) - ^l{p-^Vd,&xp{{T-T,)JBP-^Vd) 


> niVd, exp((T - Ti)JBVd) - 2, 


(2.45) 


where the last inequality is from (2.37). In the next, we will estimate gi{Vd, exp((T - Ti)JB)Vd', t € [ti,T 2 ]). 
A direct computation shows that /(£i,n) > rij2 if £i < 5 ;^, and hence B > Sri /2 := diag{\,r\l2). 
Moreover 


/i(Vd,exp((T- Ti) 7 Br/ 2 )Vrf;T e [Ti,r 2 ]) > - 1. 


V2: 


(2.46) 


71 


For s < then 


-^1 <£i,From(2.45 


2.46 1 and (2.151, we have 


p(yd,yiT,Ti)Vd\T €{Ti,T 2 \) > MVd,exp((T-Ti)7BiVrf;Te [Ti,T 2 ]) 

> /i(Vd,exp((T-Ti)7S^/2)Frf;T € [Ti,r2])-2 

^ Vnin(^) 


3. 


TT 


(2.47) 


Direct compute show that q^{Ti) - 5(£(2 V2 - £) - ((£(2 V2 - s))^ - 8 ^) 2 ) < 


8e 


3e(2V2-£) 


, then 


ln(-^) = ln(£) - )- ln(^ 2 (Ti)) > 


(2.48) 


The result is from ( |2.47p.48[ ). 

Proof of Theorem 


1.1 


Under the assumption /li(/?) = -4 - ri, from Lemma 2.6 we have For s < \s\. 


< , li{Vd,y{T,Ti)Vd\T ^ [ti,T 2 ]) > ^ln(^)-3. Similarly 

ld{Vd,y(T,T4)Vd-,T € [T4,T5]) > 4^ln(-^)-3. 

7 T ^ VI 


We have 


lu{Vd, r(T, 0)Vd-, T € [0, rj) > 2 ^ In (^) - 6. 

7 T ' 


(2.49) 


(2.50) 


The results now follows from the fact that i\{y) > fd{Vd,yVd) (for Lagrangian system, we refer §5.1 for a 
detailed discussion). This complete the proof. □ 
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3 Collision index for Planar central configurations 


This section is the main part of our paper. We give the definition of the collision index in §3.1 and we prove 
the approximation theorem; we study the basic property of collision index and we compute in detail the 
collision index on Iq in §3.2. 

3.1 Collision index 

In this subsection, we will consider the Maslov index on the half line with a hyperbolic equilibrium. This 
is similar with the case of homoclinic orbit iTTOl and heteroclinic orbit iTT^ . and a detailed study to the 
half-clinic orbits is given in lISl fT^ . 

To define fhe Maslov index of the half line, we firstly review some basic fact of heteroclinic orbits. We 
consider the Hamiltonian flow induced by 

z = JB(t)z, t€R. (3.1) 

We assume the limit is hyperbolic, meaning that 

JB{±oo) = lim JB{±t) 

t—^OO 

is hyperbolic. It follows that = Vf®V^, where Vf{V^) is the stable subspace(resp. unstable subspace) of 
the equilibria which is spanned by the generalized eigenvector of eigenvalue with negative real part (positive 
real part) of /B(±oo). Moreover, both the stable subspace T* and the unstable subspace T* are Lagrangian 
subspaces of (R^”, mo). The topology of Lagrangian Grassmannian Lag(2n) is given by the metric 

piV,W) =\\Pv-Pw II, 

where 'Pv,'Pw is the orthogonal projection to V, W and || . || is the operator norm. 

Let y(t, y) satisfy with 7 (y, y) = hn- In what follows we set y(0 := y(t, 0). Clearly y satisfies a 
semigroup property; that is, y(t, y)y(y, r) = y(t, r). For y € R, define 

K(y) - {^1^ e R^” and lim y(t, y)^ = 0), 

t—^oo 

and 

Vu{v) = {^1^ € r 2” and lim y(t, y)^ - 0). 

t—*—oo 

We remark thaf 

lim V,{v) - F; and lim F,(v) = F;. 

V—»00 V—> —OO 

It is well known that both Fi(y) and F„(y) are Lagrangian subspaces of (R^”,mo). An important property 
from im isthe following: if F transversal to Fi(0), then 

limy(t,0)F- F„F 

t—^OO 

Similarly, if F transversal to F„(0), then 

lim y(i,0)F = F;. 

t—^ — CO 

Let R- {+x > 0, X € R). We will define the Maslov index of the half line R"'' or R“. We notice 
that the discussions for heteroclinic orbit works for the half-clinic orbit. Firstly, we give the definition of 
nondegeneracy. 
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Definition 3.1. i) The linear system (113 on R is called nondegenerate if there is no bounded solution, 
ii) the linear system on R* is called nondegenerate with respect to Vo, if there is no bounded solutions on 
R* which satisfies z(0 ) € Vq- 

We observe that for the system with hyperbolic limit, all the bounded solution must decay to 0 as t —> ±oo 

LIJ. 

We firstly give the definition of Maslov index on R'*'. For, let ^ Lag{2n) and we suppose that 

the system is nondegenerate with respect to Vo> that is Vq iti ^^(0). Then 'y{t,0)Vo is a path of Lagrangian 
subspaces having limit and so, we define fhe Maslov index on R"^ wifh Fo> Vi by 

UiVuVo) := p{Vuy{t,0)Vo,t eW^). (3.2) 

In fhe case R“, recall fhaf F„(t) is a pafh of Lagrangian subspace and hm,_»_oo Vu{t) = V~, fhen for V e 
Lag{2n), we define 


L{V) ■-p{V,Vu{t),t €W). (3.3) 

We observe fhaf fhe definifion on MT does nof need fhe nondegenerafe condifion. Finally, we will define fhe 
Maslov index on R which is fully sfudied in ifT^ . Supposing fhaf fhe linear system is nondegenerafe on R, 
fhen hm,^_oo F„(f) = V~ and limf^oo F„(f) = V^. Thus we define 

i{V):=p{V,Vu{t),teW). (3.4) 

Under fhe nondegenerafe condifion, if is obvious fhaf 

i{V) - i-iV) + i+{V, U„(0)). 

We come back fo ERE. By assuming fhaf /li(/?) > we can idenfify Iq, 1+ wifh R, and idenfify 1^, 
wifh R*. Eor Vq, V safisfying fhe nondegenerafe condifions, i{Vi) on Iq or l+, /+(Ui, Uo) on or and 
i-{Vi) on /q and on Z“ are well defined, and we shall refer fo as collision index. 

Definition 3.2. The planar central configuration is called collision nondegenerate if the corresponding 
system on 1+ is nondegenerate. 

We idenfify R wifh l+, and lef Vuir) be fhe unsfable subspace. Under fhe nondegenerafe condifions, 

lim Vu(t) - V^. 

T—^±00 

Lef Vu,o be fhe unsfable subspace on l~, fhen 

lim U„,o(T) = UJ. 

r-^±oo 

Eor Vo. Li e Lag{2n), safisfying Vo iti Vi(0), fhen fhe Maslov index /+(Vo, Vi) and /-(Vi) are well 
defined. As e —> 1, we have fhe nexf approximafion fheorem which plays a key role in our paper. 

Theorem 3.3. Assuming Ti(/?) > -^, we have: (i) If V~ iti Vi, the system is nondegenerate with respect to 
Vo on Zq, and nondegenerate with respect to Vi on If then, for \ - e small enough, Vi iti jeiT' /2)Vo and 

MVi,y,(T)Vo,T€ [0,7'/2]) - Z+(Vi,Vo;Zo) + /-(Vi;Z;). 
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(3.5) 


(ii) IfV*(h V\, the system is nondegenerate with respect to Vq on and nondegenerate with respect to Vi 
on Zq, then, for I - e small enough, Vi iti 7 e(T')yJ^(T' /2)Vq and 

p{Vuye{T)re\T/2)Vo,Te[T/2,T]) = U{VuVo-,0 + L{Vul^). (3.6) 

(in) IfV^ iti Vi, the system is collision nondegenerate, and nondegenerate with respect to Vq, Vi on Zq, Zq 
correspondingly, then, for \ - e small enough, Vi itl feiT'Wo cmd 

p(Vuye{T)Vo,T e [0,r]) - Z+(Vi, Vo'Jo) + Z-(Vi;Z+) + KVuh). (3.7) 

The proof is based on a series lemmas, we firstly give the next lemma which is from ||29]| . 

Lemma 3.4. Let us consider the linear system 


X = Dx + C(t)x, (3.8) 

where D is k x k diagonal matrix and C{t) is a continuous matrix in t € [0, f], such that ||C(5)||rZ5 < e, 
for some constant e which satisfies < 1 and let y(T) be the fundamental solution of {3.8). Then for 
t€[0,f], we have 


y(Z) = (/ + O(0)exp(D0(/ + 5), 


(3.9) 


where ||O(0ll < ||.S|| < 


Proof. From lemma 6 of |[29l and for e < 1 /4, let A be an eigenvalue of D and V be the corresponding 
eigenvector.Then there exists a solution (fi(t) of (3.8 1 such that 


- ^11 < 


3£ 


1 -3£ 


Let ei,i = 1, - ■■ ,k be the canonical basis, Y(t) be the matrix defined by (fi, - ■ ■ ,(pk as column vectors. 
We define 0{t) Y{t)exp{-Dt) - I, then ||(9(Z)|| < for t € [0, Z]. Obviously, y(Z) = F(Z)T“^(0) = 
(7 + 0(0) exp(Dz)(/ + 0(0))-^ Let S := (1 + 0(0))“! - 7, then for ||O(0)|| < 1/2, we have 

ll>^ll ^ < 2110(0)11 < 


1 - 110 ( 0)11 


1 -3£’ 


which complete the proof. 


Assume that /li(7?) > -g, then P+ are hyperbolic. We recall that we set D+ = JB{P+) having the 


form Jk 


Ik 

±^Ik 


±^Ik 

-R 


4 f 


Ok h 


> P 2 = 


. An easy computation shows that the eigenvalues of D+ are real if /li(7?) > -|. 

, Tjt) with Ti < ■ ■ ■ < A/c. Let Pi = 
, P = PiP 2 - By a direct computation we get that 


Choose basis such that R is diagonalisable, that is 7? = diag{Ai 
sJl/S + RIk vJ\/S+RIk 
h -Ik 


P ^D-P = diagipi, .., pk, -pi,-, -Pk), 
where pj = -y/l/8 + Aj. Let p = max{pj, 77 “'; j = 1, • • • ,k}, then easy computation show that 

max{||P||,||p-'||)< 2 (l+ 77 ). (3.11) 


(3.10) 


Based on Lemma 3.4 we firstly prove the important lemma below. 
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Lemma 3.5. We assume that > -| and we let y(T, ti) be the fundamental solution of {2.13) with 

y{T\^T\) - hk- Then for s < sq, e < e^, we have the following estimate below 


y{T, n) = Pihk + A(t))D(t)(/2^ + S)p-\ T € [ri, T2], 


(3.12) 


where the matrices A(t), S satisfy ||A|| < yg, ||.S|| < ci£, for go.ci is constant and dependent on R and 


Proof Simple computation shows that JB(t) = D_ + C(t) with 

CiT) = 

Let W{t) = P~^y{T,T\)P, then 


0 ^ 


-q^h 

- qhn 


W{t) = (P-^D.P + P-^C(t)P)W{t). 


Let £ < Eo, where eo < 1/8 will be fixed later. It is obvious q^{T) < q{T)dT for 0 < ^ < 1. From 
equations (2.16|2.17[), we have qir) < 2e as well as \Q(t) + V2| < 2e. Then we have 


\\C{T)\\dT < £\1\\Q + ^11 + 2||<?|| + \\q^\\)dT < 6 e, 


(3.13) 


and from p.ll|) also that 


f 


\\P~^C{T)P\\dT < 24{l + fife. 


(3.14) 


Let go - (24(3 + 6 V^)(l + q)^) ^ and for g < 
From Lemma|3.41 we have 


go, we denote g = 24(1 + qfs. Then, e < eq implies 


< 1 . 


3 'sficE Cl 6 ’sficE 

IIAII < < ^g, II5II < < cig, 

1 - 3g 2 1 - 3g 

where ci = 24^ V^(l + r))^ only depend on R. This complete the proof. 


(3.15) 


Given two subspaces, graphs of two linear operators, it is possible to introduce a norm topology on 
the Lagrangian Grassmannian, equivalent to the gap topology, as below. More precisely, if F = F" 0 
a sequence of operators {L„)nen c L{E~,E"^) converges to L if and only if their graphs converge to the 
graph of L. Another important property is that, the image TV of a closed subspace V by an invertible linear 
operator T, is continuously depending on (T, V). From Lemma 3.5 P~^{Vu) = Vd and = ^n- For 

V e Lag{2n) with V iti F“, then 3Ly such that 


p-^V = Gr{Lv). 


We give a equivalent metric 

p{V,W) =\\Lv-Lw II . 

It is obvious that Ly- = Ot, since V~ iti Vi, so 3cri > 0 such that F iti Fi if ||Ly|| < cri. For cr > 0, we 
always denote 

Bp(F,(T) - {W € Lag{2n),p{V, W) < tr}. 
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Lemma 3.6. Let V = Gr(Ly), i.e. V - PV. Then for any V e U{(t) - {V : ||Ly|| < cr) with cr < I and if 
||A|| < cr/6 then (/ + A)t/ € U{2cr). 


Proof For A = 


Ai A 2 

A 3 A 4 


, we have (/ + A) 


LyX 


(/ + A] + A.2 Lv)x 
(A 3 + (/ + Af)Lv)x 


. Let 3 ^ = (/ + Ai + A 2 Lv)x 


and choose ||A|| small enough, then we have (/ + A)y = Gr((A 3 + (/ + Af)Lv){I + Aj + A 2 Ly) ^). Since 
< cr/ 6 , an easy computation shows that 


||(A3 + (/ + A4)Ly)(/ + Ai + A2Lvr\\ < 


l|A 2 ll + l|/ + A4lll|Ly|| ^ ||A|| + (1 + ||A||)||Ly|| 
/-||Ai||-||A 2 ||||Ly|| /-||A||-||A||||Ly|| 


< 2cr (3.16) 


Lemma 3.7. For any 0 < cr < 1, we let Sa- := min{eo, 2 ^ 1 - < ea-, V € B^{V^ , cr/d) and e < s^, then for 

every t e [ti, T 2 ], we have feir, ti)V € Bp{V~, cr). 


Proof From (|3.12|), p-^ye{T,Tx)V = {l 2 k+A{T))D(T){l 2 k+S)P-^V'whe.m V € Bp{y-,crlA), p-^V € t/(cr/4). 

Obviously D(T)U{crl2) c 


3.6 


Since ||5|| < Cis < then we have {l 2 k + S)P^^V e 17(cr/2) by Lemma 

U{(rl2), using Lemma 3.6 again, we have {Fk + A(t))D(t)(/2* + S)P~^V € U{cr), which conclude the 
proof. □ 


Let To be the fundamental solution on Zq as given in Equation (2.391 and T+(t, v) be the fundamental 
solution on Z+. Let cr < ^ min{p(y“, V\),p{V~, Fj)} small enough such that Bp{V~, 3cr) iti F^, Fy For this 
cr, 3cri > 0 such that Bp{V~,cr\) c Bp{V~,cr), and let Co-j be the number corresponding to crj in Lemma 

im 

Choose e < cr\ small enough such that 

max{p(To(To)Fo,F;,),p(F„(-T/J,F;,),p(T+(-Tp,0)Fi,F;,)) < cr. (3.17) 

From Lemma [ft] We have 

Lemma 3.8. For this fixed e, e < p(ye(T)Vo, V~) < crforr € (ti,T 2 ). 

Similarly, by the symmetric figure 2, for S 4 , S 5 , we have the following result. 

Lemma 3.9. For this fixed s, e < p{ye{x)y~^{T /2)Fi, F^) < crfor r e (tJ, rp. 

Proof of Theorem |3.3t To prove (i), we will compute the Maslov index p(Fi, yei^Wo', re [0, T"/2]) on 
the three time interval [0 ,ti], [ti,T 2 ] and \t 2 ,T /2]. Form Lemma 3.8 for e < s^. 


MV^i,r.(T)Fo,Te [0,Ti)) =p(Fi,To(r)Fo,T€ [0,T/J = Z(Fi,Fo;Zo). 


(3.18) 


Obviously 


fe(T)Fo, T € [Ti, T 2 ]) ^ 0. 


(3.19) 


Now we consider the path on [t 2 , T /2]. Please note that for e < s^, p( 7 e(T 2 )Fo, F„ ) < cr by Lemma 3.8 and 
p(T+(-Tp,0)Fi, F 7 ,) < cr by (13.nb, then ye{T 2 )Vo rti T+(-Tp,0)Fi, which implies 


T+(0,-Tp)r.(T2)Fo fh Fi. 
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Since j^iT - T 2 ) uniformly converges to 'F+(t, we have for e < small enough, 

MVi,r.(T)Vo,T€ [t2,7'/2]) =MVi,'F+(t,-t,JVo,t€ [-r/,,0]) = (3.20) 

The result of (i) is from ( 3.18| ), ( |3.19| ) and ( |3.20| ). 


The proof of (ii) is based Lemma [T9| and is totally analogous. 

To prove (iii), we compute Maslov index fi(Vi,ye{T)Vo',T € [0,7”]) on the five time intervals [0 ,ti], 
[ti , T2], [t2, T4], [t4, T5] and [T5, T"]. By assumption of collision nondegenerate, the system is nondegenerate 
on l+, that is limr ^+00 L„(t) = So for cr small enough, V € Ba-{V~ ), we have limr ^+00 'T(t, -ti^)V = V^. 
If we consider as the starting point, this means the system is nondegenerate with respect to V. By 
arguing as in step (i), we have for e < small enough. 

lu{Vi,ye{T)Vo,T e [T2,T4)) ^/l(Ti,'P+(T,-T/J(ye(T2)To),T e ^ 


and 


^l{Vl,ye{T)Vo,T € [T4,T5)) ^ 0, 


^i{Vl,UT)Va,T € {T5,T)) = i-iVulo), 


with (3.18 3.19 1 , we get the result. □ 


3.2 Some fundamental property of collision index 

. ( 4 f4 , 

We first compute the collision index on Zq. Recall that on line Iq, B - \ q ^ |. We can choose bases 

V 44 -R 

1 S 

f Given any two Irtik x 2mk square block matrices 


such that R = diag{ri , • • • , rk), and set Br = 

Ak Bk 


Q _r 
4 ' 


Mk = 


Ck Dk 


with k - 1,2, the symplectic sum of Mi and M 2 is defined by 


Ai 

0 

Bi 

0 

0 

A 2 

0 

B 2 

Cl 

0 


0 

0 

C 2 

0 

B>2 


Ml o M2 — 


It is clear that, B - By^ o ■■■ Br^^. We start by considering the two dimensional case. The linear systems 
z - JBrZ with the form 


y =- y + rx, 

4 

^ Q 

X = y + —X, 

4 


(3.21) 

(3.22) 


where z = (y, x)^. Assuming that r ^ 0, r > and by taking derivative with respect to r on both sides of 


(3.21 1 , we have 


Q Q 

y ^ + 

A 1 u2. V27 

- (T-x-tanh(—— 

M 8 2 
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) + r)y. 


(3.23) 

















where the second equality is from the fact that x = ^{y + jy) by (3.21 1 . Let 


•^4 8 2 


1 


V2t 


1 


) + r = -(1 - tanh ( )) + (r + -) > 0. 


(3.24) 


Lemma 3.10. Ifr > r 0, then i) for any t 2 > t\, there is no nontrivial solution of (3.23 1 which satisfies 
boundary condition y{t\)y{t\) > y{t2)y{t2)>' H) there is no nontrivial solution satisfying y(0)y(0) = 0 and 
y—>0, y—±oo; Hi) there is no nontrivial bounded solution on R. 


Proof Suppose y is solution of (3.231, then multiply by y and by integrating over the time interval [t\,t 2 \ 
we have 


r 


y{h)y{h) - y{t 2 )%t 2 ) +1 (y^ + /y^)^/T = 0. 


(3.25) 


The first conclusion now readily follows. By taking limit of ( 3.25| ) could get the second conclusion. The 
third conclusion follows from the fact that any bounded solution must decay exponential fast. □ 

Recall that Vd, are Lagrangian subspaces corresponding to the Dirichlet and Neumann boundary 
conditions. In the two dimensional case, let e\ = (1,0)^, e 2 - (0,1)^. Then it is obvious that Vd, Vn are 
the linear spaces spanned by e 1 ,^ 2 ■ An easy computation shows that z(0) = ei is equivalent to y(0) = 1 and 
y(0) = 0. In the same way, z(0) = 02 is equivalent to y(0) = 0 and y(0) = r. The second conclusion of 
Lemma 3.10| implies that for r t), the system is nondegenerate. However this is not true for r = 0. In fact, 
let r = 0 in equations ( 3.2ip.22 ). Then there is a nontrivial solution satisfies Z(0) = 02 and Z{t) 0. Thus 
we have 

Lemma 3.11. Supposing r > -1/8, the system on 1^ is nondegenerate with Vd, and it is nondegenerate with 
V„ if and only if r 0. 

We firstly consider the Maslov index on 1^, so we have the following result. 

Lemma 3.12. Suppose r > -1/8 and r Then 

UVd-,ll) = i-{Vn,ll) = 0. (3.26) 

Proof By the definition of Maslov index, we only need to show that there is no nontrivial solution. If not. 


we assume that Z(t) is the solution of the equations (3.21 3.22 1 . Then Z(t) 
y(T)y(T) —> 0. Let t\ -00 in (3.251, we have 

- yiWKti) + f (y^ + fy^)dT = 0 . 

k ) — 00 


0 as T 


-00 and then 


(3.27) 


Please note that y{t 2 )y{t 2 ) - -jy^iti) + rx{t 2 )y{t 2 ). Then for t 2 e R , Z{t 2 ) € Vd or implies y{t 2 )y{t 2 ) < 0, 
so we get the result. □ 

We continuous to compute the Maslov index on 1^. 

Corollary 3.13. Suppose r > -1/8 and r t 0, then 


i+i^n, Vd\ Iq) — 0 ) 
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(3.28) 


















and 


1 (-I’O)’ 




(3.29) 


0 


i/ r > 0. 


Proof. The proof follows from equation (2.291. Please note that Lemma 3.10 implies that there are no 
nontrivial solutions which satisfy Z(0) = e\,e 2 and Z{T) e for some T > 0. There is a crossing for 

at r = 0. An easy computation shows that the crossing form r(L„, L„,0) is positive for 
re (-g, 0) and negative for r > 0, which implies the results. □ 

To compute the collision index i+iVd, Vd\ Z2) and i+{Vd, Vn', If), we will use the Hdrmander index. We 


observe that, in the point (0, Q-), 7B^(-oo) = 


V2 

4 


r 

V2 

4 


and hence the eigenvalues 2+ = ±(g + r )2 with 


eigenvector ± + r) 2 ,1)^. The unstable subspace is spanned by and stable subspace 


is spanned by e_. Similarly, at {0,Q+), JBr{+oo) = 


iVl ^ 
4 ' 

1 ^ 
^ 4 


, the eigenvalues A± = ±(i + r )2 with 


eigenvector eX - unstable subspace is spanned by and stable subspace 

is spanned by et- The Hdrmander index could be computed by (2.351, (2.361, or we just choose simple 
Lagrangian paths connected A(0), A(l) and compute the difference. For reader’s convenience, we list the 
result below. 


s{Vd,Vn,Vd,V-) - 1, 

(3.30) 

s(Vd,Vn,Vn,V-) - 0, 

(3.31) 

' 1 i/re(-|,0), 

s{V„,Vd,V-,V:)^^ 

0, if r>0. 

(3.32) 

s{Vd, V-, Vn, f;) - s(Vd, V-, Vd, v;) - o. 

(3.33) 


The following Figures illustrate the Hdrmander index ( 3.30p.33 1, where y is the horizontal coordinate, x is 
the vertical coordinate and the anticlockwise rotation is positive rotation. 


S(Vj,,V^,V^,V‘J=0 

r>-1/8 


Vd 


.151 - 1 -^- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

-1 -O.B -0.6 -0.4 -0.2 0 0.2 0.4 0.6 O.B 1 



Figure 3: s{Vd,Vn,Vd,VX) and s{Vd,Vn,Vn,V-) Figured: s{Vd,V-,Vn,V;) and s(Vd,VX,Vd,V;) 
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Figure 5: s{Vn, Vd, F;, V^) for r € (-^,0). 


Figure 6 : ^(Fi, F^, F„ , F+) for r > 0. 


From Lemma 3.11 , the system is nondegenerate for r 7 ^ 0, so y{T)'Vd —> V^^ and 'Y{T)Vn —> Fj, as 
T —> + 00 . We have 


UiVd, Vd\ Iq) = i{Vn, Vd\ Iq) + s{Vd, Vn, Vd, V-), 


(3.34) 


UiVd, F„; /-) = /(F„, Vn, lo) + s{Vd, V„, Vn, F^). 
Thus from equations ( |3.30p.31 1 we have 
Corollary 3.14. Suppose r > -1/8 and r + 0, 

U{Vd,Vd-,ll) = \, 

and 

1 if r€ (-|,0), 


(3.35) 


(3.36) 


i+iVd, Vn', Iq) — ‘ 


(3.37) 


0 , if r> 0 . 


We come back to the higher dimension. We denote by <p{R) the number of total negative eigenvalues of 
R. It is obvious that, ^(R) = 0(ao) which is the Morse index of the central configuration a^. By property V 
of Maslov index, we have 


Corollary 3.15. Supposing /li(7?) > -1/8, R is nondegenerate, then we have 

UVd', ID = 0, L(F„; ID = 0, 


(3.38) 


and 


i+{Vd, Vd', Iq) = k, i+{Vd, Vn', ID = 


i+iVn, Vd', 1 q) — 0, i^{Vn, Vn', ID ~ 


(3.39) 


(3.40) 


We then consider the linear system on I+, and recall that under the collision nondegenerate condition 
/(F) is well defined. 
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Proposition 3.16. Under the collision nondegenerate condition, we have 


i{Vn,h) = i{Vd\h) + cp{R) (3.41) 

Proof. Please note that i{Vn', 1+) - i(Vd', 1+) + s{Vn, Vd, V~, Vf), we only need to show 

s(Vn,Vd,V-,V:) = f{R). (3.42) 

We choose basis such that R = diagifi, ■ ■ ■ ,rn) and consider the 2-dimension linear system z - J 2 B 2 {r)z', 
(13.421) is from (|3.32 1 . 


Now we give the proof of Theorem 1.2 


Proof. Since R is nondegenerate, y„, iti P*, from (3.7 1 , 

\im^i{Vd, 7e{T)Vd, T e [0, T]) = /.(P^; Z+) + U{Vd, P^; Zq) + Z+), 

e —>1 

from p.38p.39[ ), we have 


\imii{Vd, ye{T)Vd, T € [0,7-]) = k + i{Vd\U). 

e—>l 


Then, (1.6 1 follows from Lemma [2!4| and the fact that 7e(27r) = jeiT). 

Similar, 

lim//(P„,y,(T)P„,T € [0,r]) = UVnUl) + /+(P„, P„;Zo) + KVnJr), 

e-^i 


([HI) is from p.38| ), (| 3.40H3.4T] ) and Lemma [24 


(3.43) 


It is worth noticing that in several cases we can’t analytically compute the collision index and for this 
reason, we now introduce a numerical method which can be useful in this situation. 

We first consider the case for R"*". We choose P e Lag{ln), such that B{t)\v > 0 for t € R*. Then the 
crossing form r(A(T), P, r) > 0 and we have 

piV,y(.T)Vo) = ^ v ( t ^), 

0<Tj<OO 

where v{Tj) = dim P n yirflV^. For the Lagrangian system, we can always choose Vd - V and then 
= 7 h > 0. We can get the Maslov index from the Hormander index, in fact 

/i(Pi,r(T)Po,Te [0,r]) =/i(Prf,7(T)Po,[0,r]) + 5(Pi,Pd,Po,r(r)Po). 

Under the nondegenerate conditions, lim 7 ’_»oo 7(7 )Pq = Vf, then we have 

Ai(Pi,7(T)Po,T e R^) - KVd,y{T)Vo,T e R+) + ^(Pi, Vd, Po, V^). 

The cases of R“ and R are similar, so we just make -7 be the starting time, where 7 > 0 is large enough. 

Remark 3.17. For computing the Maslov index p{Vd, 7(t)Po), we start by choosing a basis [fi (0), • • • , fn{0)} 
of Po and, by using a numerical integrator, we get fk{t) '■= y(t)fk{0). For j = I, - ■■ ,n, let ej the basis of 
Vd, t i-> M{t) be the path of In x 2n matrices defined by M(t) := (ei,--- ,e„,fi{t),--- and we set 

f(t) = det(M(t)). Then p{Vd,y{t)Vf) is equal to the total number of zeros of f{t). Since t |/(t)| is 
exponentially increasing, this method works well for not so large time. Instead, we use the robust numerical 
algorithm based on exterior algebra representation. We refer the interested reader to the following papers 
(013 ID/ in which the authors compute the Maslov index for homoclinic orbits. (Cfr. Section §6). 


24 














4 Collision index for brake symmetry Central configurations 

In the next, we will consider the case that the central configuration with brake symmetry. We start with the 
following Definition. 

Definition 4.1, The central configuration with normalized Hessian R is said with brake symmetry if there 
exists a kxk symmetry matrix N which satisfies = Ik, NJ = -JA^, RN = NR. 

To our knowledge, the Lagrangian configuration ifTTl , Euler central configuration ETl and the \ + n 
central configurations have the brake symmetry property 041 . HOl . It could be interesting to find an example 
wifhouf having fhis symmefry property. 

In fhe brake symmefry case, lef N - diag{N, -N), and denofe 

g : x{t) Nx{T - t). (4.1) 

Obviously, g^ = id and g ■ -4^ = “4^ • g. From fhe facf of q{T) = q{T - r), Q{t) = -Q{T - r), easy 
compufafions show fhaf -t) - B{t)N, and consequenfly 

gB - Bg. (4.2) 

Lef E- = ker(g + I), fhen 

ker(-7|- -B) = ker((-4£ - B)|^J © ker((-4|: - B)|^_). (4.3) 

Moreover, by fhe generalized Boff-fype iferafion formula for Maslov index lIT^ f Thl.l) or |[26l . we have 

/i(r) + k = iu{V^{N),y{T)V^{N),T e [0,T/2]) + fi{V-{N),y{T)V-{N),T e [0,r/2]), (4.4) 


/-i(r) - ti{V^my{T)V-{N),T e [0,r/2]) +/i(I/-(iV),r(T)E+(fV),T € [0,7'/2]), (4.5) 

where V-{N) - ker(fV + hn), k - 2n-4. Similarly, we can decompose of Dirichlef and Neumann boundary 
condition as follows 

BiVd, y{T)Vd, T € [0, rj) = y{T)Vd, T € [0, r/2]) + id{V-{N), y{T)Vd), t € [0, r/2]), (4.6) 


B{V„,y{T)Vn, T € [0, T]) = y(T)Vn), T € [0, r/2]) + iu(V-(N), y{T)Vn), T € [0, r/2]). (4.7) 


For fhe iferafion formula of brake orbifs we refer fhe inferesfed reader fo 11231 . Then we consider fhe collision 
orbif, and lef TC be fhe space of bounded solution of z = JB{t)z on /+, and TC© be fhe space of bounded 
solution on which satisfies z(0) € V-{N). Similar fo (4.31, we have 


Lemma 4.2. For the brake symmetry central configurations, on l+, we have 


7C = 7C+©7C-. (4.8) 

A. A. A. A. A. 

Proof. Please note fhaf on Z+, NB{-t) = B{t)N, if z(r) is one solution fhen Nz{-t) is anofher solution. Lef 
z+(t) = |(z(t) ± ZVz(-t)), fhen z+ € 'K^, which implies fhe resulf. □ 
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Obviously, there is standard brake symmetry on Iq, that is Nq - diagih, -4)^ then V‘(iV^o) - and 
V~{No) - Vn- Let TCo be the space of bounded solution of z = JB{t)z on Zq, and 1^+ be the space of 
bounded solution on which satisfies z(0) € V-{No). We have on Zq, 

Wo = ’Xt©’X‘°. (4.9) 


On Zq, the nondegenerate condition is clear. Now, from (4.91 and Lemma [3.1 1[ we have 
Proposition 4.3. The system on Zq is nondegenerate if and only ifR is nondegenerate. 

Since R satisfies brake symmetry, let V- - ker(At + 4), so it is obvious dimV- - |. Let V- - JV- © V- 
be symplectic subspace thus = L'*' © V~. Set and denote B+ - 


k 

2- 

^ Ik 
V 4 2 


, we get 
-R^ ^ 


4 = B+ o B_, and the fundamental solution satisfies 'Pq = 'Polf/+ o 'Loly-- Furthermore we have 


i{V, W; Zq) - mv,,W\vT,lo) + Z(F|w, W|y_; Zq), 
for V,W = V-{N) or Vd, V„. Please note that 

= Vd, = V„, F-(4)|y. - W(4)|y- = 


(4.10) 


(4.11) 


From (3.39 1 -(3.40 1 , we have 


U(V\N), Vd-, Zo) + h(V-(N), Vd-, Zo) = k, 
U(V^(N), Vn-, Zq) + Z+(V-(4), Z") = 2f{R), 


i^iV^m V\N)-, Zo) + Z+(V-(4), V-{N)-, lo)^k + fiR). 
U{V-{N), V^{N)-, Zo) + UiV^iN), V-{N)-, Z^) - f{R). 


(4.12) 

(4.13) 

(4.14) 

(4.15) 


For the brake symmetry central configurations, we get the following approximation theorem. 

Theorem 4.4. Let A\{R) > - ^ be nondegenerate with brake symmetry property, and satisfying the collision 
nondegenerate conditions. Thus we have 

limp{Vd,7e{T)Vd,T e [0,T]) = k + Z_(F44);Z;) + UV\N)-,lf), (4.16) 


limMV„,r,(T)F„,T € [0,r]) = 2<P{R) + L{y-{N)-,l-) + Z_(F^(ZV);Z;), (4.17) 

lim Lfje) = lim ifje) = cf>{R) + i-{V-{N)-, If) + i-{V\N)-, If). (4.18) 

e—^\ 4— 
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Proof. Since R is collision nondegenerate, from (4.8 1 , the system is nondegenerate with respect to V-{N) 
on l~, and also from the condition that R is nondegenerate, then the system is nondegenerate with respect to 
V-{N), V„, Vd on l~. From (i) of Theorem : 


3.3 


we have 


lim/i(F,f,(T)lT,Te [0,r/2]) - /(F, IT; Z^) +/(V;/;), 

e-»l 


for V, W is V-{N), Vn, Vd- From (4.4 4.7 1 , we have 


lim^i{Vd, 7 e{T)Vd,T e [0,r]) = P{V^{N), Vd', Zq) + P(V-(N), Vd; Zq) + Z_(F+(ZV); Z;) + Z_(F-(ZV); Z;), 

e-*l 

limid{Vn,7e(T)Vn,Te[0,T]) = U{V^{N),Vn;lo) + U{V-{N),Vn;lo) + L{V*{N);l-) + L{V-{N);l-), 

e-*l 

lim Zi(r,) = Z+(F+(ZV), F+(ZV); Zq) + L(V^{N); Z;) + Z+(y-(ZV), V-{N); Z") + L(V-(N); Z;) - k, 

e —>1 

lim Z_i(f,) - U{V-{N), V^{N); Z^) + Z_(F+(ZV); Z;) + U{V^{N), V-{N); Z") + Z_(y-(ZV); Z;). 

e —>1 

Then ( 4.16||4.18 1 is from ( |4.4H4.7| ) and ( |4.12H4. 15] ). 

Comparing of (1.6 1 in Theorem |1.2| and (|4.16|), we have 


Corollary 4.5. Let AfR) > -1/8. If R is nondegenerate, has the brake symmetry property and if the collision 
nondegenerate condition is fulfilled, then we have 


Z_(y-(ZV); If) + L{V\N); If) = i{Vd; Z+). 


(4.19) 


It is clear that i(Vd; If) ^ 0. We shall now prove that i-{V-{N); f) is also nonnegative. We consider the 
Maslov index on R“, supposing the system is nondegenerate with respect to Vi, and Vi iti V~, then for -tq 
large enough, 

p(.Vu Vu{t),t e (- 00 , 0 ]) ^ /l(yi,7(T,To)y„(To),T € [to,0]). 

From the property (III), (IV), for 7 (- 5 ,0), 5 € [0, oo) 

MVi,r(T,To)y„(To),T e (to,0)) =p{7{-s,0)Vi,V~,se [0,oo)), 

we have 

Z^(Vi,y„(r),T € (-00,01) = -p(.V~,7(-T)VuT e [0,oo)). 

By the nondegenerate condition, we have lim 7 '_»oo 7{-T)V\ - Vs, then we have 

p(Vu y„(T),T € (-cx.,0]) = s(Vd, V-, Fi, v;) -p{Vd, 7 i-T)VuT e [0, +CX,)). 

In the case of ERE, an easy computation shows that 


—'F+(-t) = -7B(-t)T'+(-t), 
dr 


(4.20) 
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where 'P+(t) = 'F+(t, 0) is the fundamental solution on /+. If the central configuration satisfies the brake 
symmetry, that is - B{t)N, then, direct computation shows that 

'F+C-t) - iV'P+(T)iV, T € [0, oo). 


So we have 


IJiVd,^^(-T)Vi)=luiVd,N^'^iT)NVuT€[0,c^)). 

Please note that if -B{t)\v < 0 for t € R’*', then the crossing form r(A(t), Vd, t) < 0, we have 

0<ry<oo 

where v(t;) = dimP^ n 


Please note that, in the case Pi = vi 0 where vl € Vd, € P„ , from (13.331), we have 


s{Vd,VZ,Vi,V:) = Q. 


Since V-{N) is a direct sum of Dirichlet Lagrangian subspace and Neumann Lagrangian subspace, by (4.11 1 , 
we have 


Lemma 4.6. On L, we have 


i-{V^{N))= Y 


(4.21) 


where v{tj) - dim Vd n 'i'+(-Tj)V-{N). 


5 Applications 

We give applications for the collision index. In subsection §5.1 we consider the ERE of minimal central 
configurations and prove some hyperbolicity results. At §5.2, we study the stability of Euler orbits. 

5.1 Minimal central configurations 

In order to give a hyperbolic criteria, we first review some results on Morse index. Consider the linear Sturm 
systems 

- ^(P(t)y + Q{t)y) + Q^{t)y + R{t)y - 0, (5.1) 

dt 

as P,R, Q are continuous path of matrices in R^” and satisfy P{t) > 0, R{t) = R{t)^. This linear Sturm system 
( |5.1| ) corresponds to the linear Hamiltonian system 

z = 7R(0z,zeR2”, (5.2) 


where 


' p-\t) -P-\t)Q{t) 

.-Q^P-\t) Q\t)P-\t)Q{t)-R{t) 
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B{t) = 


(5.3) 





Let 


1 r 

L{t, x(t), x{t)) = + Qx) ■ X + Q X ■ X + Rx ■ x), (5.4) 

and Tix) - {L{t, x(t), x{t))}dt on W^’^([0, T], C”). We denote 

D{oj, T) = {xe W'’2([0, r], C”), x(0) - cox{T)}, coeU. 


Obviously, 

Wo’^([0, r],R") c D{ 0 J, T) c w'’2([0, r],C”). 

Let £ = £'"(0) := -£(R(0£ + 2(0) + 2^(0ff + ^(0> and more precisely, set £n,£oj,£d to be the 
operator with form £ under the Neumann, m and Dirichlet boundary conditions separately. Let Ak{£) be the 
^-th eigenvalue of £. From the monotonicity property of the eigenvalues lIT^ . we have 

AkUn) < AkUc.) < AkUd). (5.5) 


Let (j) be the Morse index of £ which is defined to be the total number of negative eigenvalues, which is 
equal to the dimension of maximum negative definite subspace of T. Let (l>d, (pa, (pn be the Morse index of 
£d, £w, £n separately. From ( |5.5[ ), we have 

(pd ^ ^ (pn- 

Proposition 5.1. The system is hyperbolic if (pn = fd and £n is nondegenerate. 

Proof. Please note that the system is hyperbolic, that is, cr{y{T)) n U = 0 is equivalent to £a which is 
nondegenerate for Va» € U. Supposing ko - (p„ = (pd, we have AkQ{£a) < 0 by (5.51. On the other hand, £„ 
is nondegenerate which implies AkQJr\{£n) > 0 and hence AkQ+ii£cj) > 0, which implies the result. □ 


Please note that = 0 implies (pd = 0, so we have 
Corollary 5.2. The system is hyperbolic if £n > 0 . 

From Theorem 1.2 of |[T8]| or P172 of |[24l . we list the relation of Morse index and Maslov index below. 
Lemma 5.3. Let y be the fundamental solution 0/(0, then we have 

(pM - Liy), Va{£) = Vaiy), dco € U, (5.6) 


(Pd{£)+n^p{Vd,yVd), (pn{£) = li{Vn,yVn). (5.7) 

Proof of Theorem |1.3| Please note that the central configuration ao is non degenerate minimizer which 
implies /li(/?) > 0, i.e. (p(R) = 0 and R is nonsingular. Under the collision nondegenerate condition, from 
(|L6|)-([L7|), for 1 - e small enough 


p{Vd,ye{t)Vd,t e [0,2n-]) - k = \mip{Vn,ye{t)Vn,t € [0,2n-]). 

e-il 


(5.8) 


From ( |5.7[ ), we have fd = fn- The nondegenerate of £„ is from Theorem 1.2 so the result is from Proposi¬ 
tion [SH] □ 

A typical example is the Lagrangian equilateral triangle central configuration. It is obvious that R = 
diag{{3 ■+■ yj9 -/3)/2, (3 - yj9 -/3)I2) satisfies the brake symmetry with N - diag{\, -1). This fact had 
been used to decompose the -1-degenerate curves in ITTl . It is proved in ifTTll that for any € (0,9], I - e 
small enough, £n is positive, and consequently hyperbolic. By the approximation formula ( |4.17[ ) and the 
nonnegative property (|4.21|), we have 
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Proposition 5.4. If the Lagrangian central configurations is collision nondegenerate, then i-{V (N)', If) - 
i-{y^{N)\ If) = 0, and hence i{Vd', If) = 0. 

We continue by studying the case of strong minimizer, so please note that a central configuration is 
strong minimizer if it satisfies /li(/?) > 1. The nexf lemma is imporfanf in fhe proof of Theorem |1.41 

Lemma 5.5. (See Proposition 2.) If 6 > I,to e U, then fi{e, d) = -1 + is positive operator 

for all e € [0, 1) on its domain Dfcji),ln), where Dn(to,2n) = {y e W^’^([0,27r],Cf\y(27r) = ojy(0),y(27T) = 
ojyiO)}. 

Now we can proof Theorem |1.41 
Proof For fhe ERE, we have 


d^ ^ d R 

-C = ~ -TT- 

dr- dt I + e cos(t) 


fhen £> X, ■- 


-^h - 2ikiipf + i+‘/fos(f) • need fo show X > 0 wifh domain Diftjj,2n) for any 


CO eU. Eef Rlt) = 


cos(t)4 -sin(04 ] 
sin(t)4 cos(t)4 / 


Rln^ 


dt2 k ^ I ^ g COS(t) ’ 


Since Ti(/?) > l,we gel fhe resull from Eemma[53| 


□ 


If is obvious lhal for fhe slrong minimizer if if is collision nondegenerale, fhe approximation fheorem 
implies i(Vd, If = 0. 

In fhe special case, fhe ERE of Eagrangian cenlral configuralions is hyperbolic for 4 > 8, e € [0,1), 
which had proved direclly in lIMl . As anofher example, we consider fhe 1+3 cenlral configuralions, and lei 
m\ - m 2 = m^, - I and mo = me, Ihe essential pari R - 1a,+ 1) wilh 


1 

2 

0 

3 V3mc(3+mc) 

2(1+V3m£:) 

0 


0 

1 

2 

0 

3 fmcO+nic) 

2 ( 1 +V 3 mr) 


3 d3u(3+mc) 
2{l+fnic) 

0 

V3(3+m,) 

2 ( 1 +V 3 mc) 

0 


3 d3>nc(3+mc) 
2(1+V3m<;) 

0 

V 3 ( 3 +mc) 

2{1+V3m<;) 


(5.9) 


was computed in IISTl . Please note lhal Ihere is a lypo in (39) of OTII . Eel 


]_ — 3 V3mc(3+mc) ' 

<p) _ 2_ + 2(1+V3m,) 

^ _ 3 V3mc(3+mc) fO+nic) ’ 

, 2(1+VSm^) 2(1+VSm^) 

Ihen D = ID- o 2)+. Obviously Ihe eigenvalues A± of 2):; is same, direcl compulation shows lhal 


A±(mf) = -(1 + V3m, 


,) '[V3mc 


+ 


3V3 + 1 


(27(m: 




(5.10) 
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Obviously A+{mc) > 0 for rric € [0, +oo). Let mdO) = md-l) = ^^ 249 ^ ’ A-imdO)) = 0 and 
A-inici-l)) = -1, moreover 

A-(mc) >0 if nice [0, mdO)), 

■ -1 < A-(mc) <0 if nice (mc(0), mc(-l)), (5.11) 

< A-{mc) < -1, if nic e {md-l), + 00 ). 

Since = 1 + A-{mc), it is obvious that Theorem |1.3| and Theorem 1 1.3 1 imply Corollary |1.51 

Remark 5.6. Inspired by the recent results obtained in M 71/ . we conjecture that the nondegenerate minimal 
central configuration is collision nondegenerate and satisfies 

i{Vd\h) = 0. (5.12) 

In a private communication with the first name author, prof Y. Long posed the following conjecture. 

A smooth T-periodic non-collision solution of the planar N-body problem, with N > 3, is a 
smooth global minimizer of the action functional on the space of all T-periodic orbits having 
non-trivial winding number if and only if it is an elliptic motion corresponding to the global 
minimal central configuration of the potential restricted to the inertia ellipsoid. 

We observe that. Long’s conjecture implies that, for any e e [0,1) the Morse index for the ERE of any non¬ 
degenerate minimal central configuration is 0. In the case of brake symmetric central configurations. Long’s 
conjecture, imply our conjecture in the collision nondegenerate case. This conjecture is still open, and in the 
case e = 0, an interesting result from Chenciner and Desolneux shows that the minima of the action on the 
zero mean loop space, is the relative equilibrium corresponding to a minimal central configuration. UIV . 


5.2 Stability analysis of Euler orbits 


The Euler orbits have been studied in ESi . ETl . in this case, R = diag(-5,2S + 3), where 6 e [0,7] only 
depends on mass mi,m 2 ,m^. Please refer to Appendix A of Il28l for the details. Although there is no 
physical meaning for d > 7, we will assume d > 0 to make the mathematical theory complete. 

We will use the index theory to study the stability problem. Let yg^e be the fundamental solutions of 


Sf) which is given by (2.2), that is yg^e = J 2 lB{t)yg^e, t £ [0,2;:], yg^eiO) = I 4 . For d = 0, the system 
degenerates to the Kepler problem, and it has been studied in |[T9ll . which proved that yo,e(2n') with normal 

form I ^ ^ I o /2 and 

01' 


Liyo.e) = ■ 


0 , 

2 , 


if or = 1, 
i/m€U\{l}. 


(5.13) 


The Maslov-type index in the case d > 0 has been studied by Long and Zhou liT/ll . then we review firstly 
their results. For any j e N, there exists 1-degenerate curves Ty = Gr{(pj{e)), and we also let Fq = Gr{(po{e)) 
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with (fiQ{e) = 0. Then only degenerates at iTy and dimkeT{yse{2n) - I 4 ) - 2 if {5,e) € U“ ^Ty The 
Maslov-type index satisfies 


iii76,e) ^ 2j + 3, if (fife) <6 < (pj+i{e), y € N U {0). 


(5.14) 


Similarly, for Vy € N, there exists pair -1-degenerate curves T| = Gr{ij/'j(e)). 

Let = min{il/*j{e),4f~{e)] and = max{il/*j{e),tjj~{e)}. Moreover, we set iJ/q = iI/q = 0, then for 
e N we have 


i-i{7s,e) 



if Se 


(5.15) 


27 + 1, i/de(^;,«A']. 

Direct computation shows that i/^|(0) = but it is not clear if, for e > 0, there exist other intersection 

points. There is a monotonicity property for Maslov-type index, that is for m e U 


L{76i,e) < icoi7d2,e), if < ^2. 


(5.16) 


The ± 1 degenerate curves satisfies 

0 < < 1 A 2 ^ 1 A 2 < • • • ^ < ‘fij < 1 A 7 ^ <A;- < • • • ■ (5.17) 

Moreover for fhe region befween fhe ±l-degenerafe curves, ys^ei'^n) is ellipfic-hyperbolic and for fhe region 
befween fhe pairs of -1-degenerafe curves ys,e{2n) is hyperbolic. 

As a confinuous work of Long and Zhou ll27l . we use fhe near-collision index fo sfudy fhe limif case. 
From Theorem 1.1, we have 

Theorems.?. If 5 e (1/8,7], let e = ^minl^^, 1/8), e = Fore < e^, we have 

/i(y.)>-(d-^)hn(^)- 6 . (5.18) 

TT 8 ^ 

Hence 


sup{lime^i(7!7(e), limg 




< 1 / 8 . 


(5.19) 


Proof Since/? = d/ag(-d,25-i-3), fhenfor d > 1/8,/li(/?) < -1/8. Now equafion (5.181 follows from ( [T3] ) 
of Theorem 1.1 To prove (5.19), lef d/ = limg^upfe), fhen fhere exisfs ei ^ \ such fhaf lim/^oo ^fei) = 6j. 


If 6 j > 1 /8, fhen we choose dg € (1 /8, Sj), for I large enough, dg < (ffei), by fhe monotone properly (5.161, 
so we have iiiys^^ef <2j+l, which conlradicl to (5.181. The proof for lime_»ii/^*(e) < 1/8 is similar. 


Lef N = diag(\,-l), fhen NR = RN, we will compute fhe collision index for d e (0,1/8) by fhe 
decomposilion properly. By fhe brake symmelry, from (|4.3[) we have 


d/mker(y5,^(27r) -Hi) - dim{V~{N) n ys,e{2n)V^ (N)) + dim(V*(N) n y5,^(2;r)F“(/V)). 


We always sel lo be fhe degenerale curve in fhe sense lhal V^{N)C\yg^eQ-Ti)V (N) nontrivial and similarly 
to be the degenerate curve in the sense that V~{N) D y 5 ^e{ 2 n)V*{N) nontrivial. 
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We get the collision index on numerically from (4.21 1 and the step in Remark 3.17 With the help of 
matlab, we have 

Numerical result A: For the Euler orbits is collision nondegenerate for S e (0, 4), and on l~ 


L(V^(N);0 = L(V-(N);i;)^l. 


It is obvious that 0(7?) = 1 for d € (0, 4), then from (4.18 1 , we have 


Corollary 5.8. For 5 € (0, 4), under the condition of numerical result A, we have 


lim/i( 7 ) = lim/_i( 7 ) = 3. 

e-^\ e-A{ 


From (3.39 


3.401, easy computation shows that for 6 e (0, g), 

i{V\Nl v-(Ny, / q ) = 1 , i(V-{N), V^iNy, /+) = 0 . 


So for 1 - e small enough, 


iaiV^iNyys,eV-{N)) = 2, 


(5.20) 


(5.21) 


(5.22) 


KV-(Ny7s,eVHN)) = 1 . 

Theorem 5.9. Under the assumption of numerical fact A, for the ±l-degenerate curve, we have 

\imipj{e) ^ limi/^7 j(e) = 1/8, for j e N, 

e—^\ e—^\ 

\mnp\{e) = 0, \mif1{e) = 1/8. 


Proof To prove (5.24i, from Theorem|5.7[ we only need to show 


inf{lim (file), lim tJ/~j^fe), j e N) > 1/8. 

e—^ 1 e—^ 1 


(5.23) 


(5.24) 


(5.25) 


(5.26) 


The proof of (5.26l is similar to (5.191. Let Sj = ]\m(fij{e). If dy < 1/8, then we choose ei —> 1, such that 


(ffei) — > dy. Choose e < 1/8 - dy, for I large enough, ipfei) < 6j + e, so we have 

h{76j+s,e,) > - 27 + 3 , 


(5.27) 


which is contradict to (|5.21|. It is totally similar that lm(/^*^j(e) > 1/8 for 7 e N. 


e—*l 


Direct computation shows that for d € (0,1/8), p{V iN), 7 s,oV*{N)) = p{V*{N),ysflV (N)) = 1. By 
monotone property. 


p{V\Nyys,eV-{N)) = 


1, i/d€(0,«A4], 

i 2, i/d€(<A|,^J]. 


(5.28) 


From (5.22 1 , we get limg^i i/^|(e) = 0. The proof for lime_>i i/^j (e) = 1/8 is from (5.231 and the step is 
similar. □ 
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This theorem shows that the system is hyperbolic for 6 e (0, g), and \ - e small enough. To explain the 
results, we use the following pictures which are taken from If29l . 



Figure 7: Stability bifurcation diagram. Figure 8: A magnification of Figure 7 for 1 - e small. 


6 Numerical results for collision index 


As shown in Remark 3.17 in order compute the collision index, we only need to count the zeros of a 
determinant function. We use the exterior algebra representation from 171013 to do the computation. For 
reader’s convenience, we give a brief review in the four-dimensional case, here. 

Consider the linear system 


X - A(t)x, X e T e [0, -i-oo). 


(6.1) 


where A(-i-oo) is hyperbolic. Let A^(R^) be the vector space of 2-vector space in R"*. Supposing ej,] = 
1,... ,4 is basis of R^, then ei = ei A 62 , 62 = e\ A e^,, = ei A 64 , 64 = 62 A e^, 65 = 62 A 64 , e(, = A 64 
is basis of A^(R^). There is a induced system from (|6.l| 


y = A(2)(T)y, y € A^CR^). 


( 6 . 2 ) 


Suppose A = (atj), then A*^^^ could be expressed by (aij). (Cfr. 0 Equation (2.8)] for the expression). Let 
cr be the sum of the eigenvalues of A(cx3) with positive real part. Let y(T) = e~°’’’y(T), then 


f=(AAr> 

CIT 


cr/4)y. 


(6.3) 


To compute the Maslov index fi{Vd,yiT)Vo), we choose a basis ^i(0),^2(0) of Vq, and let y(0) = y(0) = 
^i(O) A ^ 2 ( 0 ) = Z/=i yA^)^i- Then %t) could be computed by matlab from Equation (6.3 1. Let y be the 


fundamental solution of ( |6.1| ), then y{T)Vo could be expressed by y(T). We choose ei,e 2 to be the basis of 
Vd, then it is obvious that Vd n y{T)VQ is nontrivial if and only if ei A ^2 A y(T) = 0, which is equivalent to 
y(,{T) - 0. So we can draw the picture of y(,{T) and count the number of zero points to get the Maslov index. 
We will compute i-{V-{N); Z“) for Euler and Lagrangian orbits. Lrom Lemma 4.6 we only need to 


count the points of Vd n y{-T-)V-{N). Lrom (4.20|, the linear system with form with x(t) = -JB{-t)x{t), 


let A(t) = then we can get A^^^(t). We choose e\,e 4 to be the basis of V*{N) and 62 , to be the 
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basis of V~{N). Let be the solution of (6.31 with initial condition 3)(0) = and y~{T) be the solution 
with initial condition j)(0) = € 4 , then L(L-(iV); l~) just is the zero points of Jg (r). 

We firstly give some numerical pictures for Euler orbits: 


0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 
0 

- 0.02 

0 5 10 15 20 0 5 10 15 20 




Figure 9: ygC^r) for d = 0.1. Figure 10: yg(T) for d = 0.1. 

It is obvious that there is only one zero point in Figure and Figure and we have computed it for 
many value of d e (0,1/8) and for time large as r = 1000. All the pictures shows that there is only one zero 
point. This is why we gave Numerical result A. 



0.16 

0.14 

0.12 

0.1 

0.08 

0.06 

0.04 



- 0 . 04 '-^-'-'-'-' 

0 20 40 60 80 100 


Figure 11: y^Cr) for d = 1. Figure 12: y^ir) for d = 1. 

If d > 1/8, Theorem 1 1.1 1 shows that the collision index is infinity, which is corresponding to the picture 
of Figure [TT] and Figure [^that the number of zero points growth in direct proportion to the time. 

For the Fagrangian orbits. Proposition |5.4| shows that the collision index is zero, which corresponds to 
the following pictures showing no zero point. 
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Figure 13: Jg(T) for {3 = 6. Figure 14: y^ij) for f3 = 6. 

It is quite difficult in a concrete situation to establish if an orbit is collision non-degenerate. However if 
the the collision index depends on one parameter having a jump, then there is a collision degenerate point. 
During our computations, we found that the Kepler case (6 = 0) is collision degenerate. We guess that every 
nondegenerate central configuration satisfying the condition /li(7?) > -1/8 is collision nondegenerate. 
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